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Abstract 

We apply a simple analytical method based on a unitary transformation to calculate the Bloch- 
Siegert (BS) shift over the entire driving-strength range. In quantitative comparison with the 
numerically exact BS shift obtained by Floquet formalism as well as the previous BS results, we 
confirm that our calculated results are not only accurate in the weak-driving regime but also correct 
in strong-driving limit. In the intermediate strong-driving regime, the calculated values of the BS 
shift are nearly the same as the exact ones. It turns out that our calculation for the BS shift is beyond 
perturbation. Meanwhile, we demonstrate the signatures caused by the BS shift by monitoring the 
excited-state population and the probe-pump spectrum under the experiment accessible conditions. 
In particular, we find that when the driving frequency is fixed at the transition frequency of the 
system, the lineshape of the probe-pump spectrum becomes asymmetric with the increase of the 
driving strength, which may be verified experimentally. 
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I. INTRODUCTION 


The physics of driven quantum system^ as an attractive topic in quantum physics, has 
been widely studied for several decades l|-l6j. At present, the study of such systems is 


atoms such as superconducting circuits 


13 


significant for quantum information processi^ j?- 12| and renewed in the context of artificial 


16| . The prototype of driven quantum systems 


is the Rabi model describing a two-level system (TLS) driven by a harmonic driving (we set 

h=i). 


1 A 

H{t) = -cuoo-^ + — cos(a;t)aa, 

1 A 

— A (J_ -f e 


(^+) 
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(e cr_ -I- e cr+j. 
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where <Jx,y,z is the usual Pauli matrix and a± = {ax ± iay)/2. coq is the transition frequency 
of the TLS. A and u are the amplitude and frequency of the harmonic driving, respec¬ 
tively. Usually, one invokes the rotating-wave approximation (RWA), i.e. the omission of 
the counter-rotating (CR) terms -1- 6*^^^+) 17j. The RWA is valid in the weak 

driving and resonance cases, but breaks down in the strong driving regimes [2]. It has been 
recognized that the CR terms lead to the shift of resonance position and additional beats 


in the time evolution 


QBi. 


Particularly, the shift of resonance position from resonance 
frequency of the RWA case is called the Bloch-Siegert (BS) shift: (Jcubs {l, 2|- We focus on 
the BS shift in the present work. 

The BS shift has been considered extensively in the cases of both classical and quantized 


fields 


flis, 


18 


2n|. For the classical field, it is convenient to give the BS shift by the Floquet 
formalism developed by Shirley j2|. For the quantized field, a fully quantum mechanical 
description was provided in the so-called dressed-atom model that combines the atom and 
driving field 20|. All in all, most of pervious works, which are based on perturbation 




Hence, there is not yet a 


calculation, give the well-known BS shift in 6th order of A 
simple analytical method to produce the BS shift in the entire range of the driving strength. 

Apart from the theoretical calculation of the BS shift, the effect of the BS shift has 
been studied before, such as an asymmetric Antler-Townes profile in a driven three-level 

application to 


system 


2l| . the shift of the sidebands of resonance fluorescence 


in 

a c 

22, 

23 


the astrophysical determination of the fundamental constants variation 


241, etc. More 
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recently, we showed that the BS-type correction plays an important role in the relaxation 
and dephasing processes j^. In addition, the BS shift has also been experimentally observed 
in strongly driven artihcial atoms {l5, 16|. 

In this paper, we calculate the BS shift over the entire driving-strength regime (0 < 
AjijjQ < oo) by a simple analytical method based on a unitary transformation The 
approach has been used to study the driven tunneling dynamics of the TLS in our previous 
work It turned out that our method provides almost the same description as the numeri¬ 
cally exact treatment over a wide parameter regime of interest. Moreover, we prove that the 
BS shift up to the 4th order in A given by our method is the same as that obtained by the 
Floquet approach Q. In this work, we show that both small and large BS shifts can be uni¬ 
formly and correctly evaluated from the derivative of effective Rabi frequency we derived, 
in comparison with numerically exact results and previous analytical results. Moreover, 
we illustrate that the signatures of the BS shift can be monitored in the emission process 
and probe-pump spectrum under the experiment accessible conditions. When u = uq, we 
show that the sidebands of the non-RWA probe-pump spectrum are generally asymmetric. 
The asymmetry results from the BS shift and can be enhanced with increasing the driving 
strength. On the other hand, the non-RWA probe-pump spectra become symmetric only 
when oj = uJo + Subs- According to these properties, one may check the signatures of the BS 
shift experimentally. 

II. UNITARY TRANSFORMATION 


The wave function |'k(t)) of the TLS satisfies Schrodinger equation 

d 


Hit) - i 


dt 


!»(«)) = 0 . 


We transform the above equation as follows: 

d' 


A(t) 


H{t) - i 


dt 


[e-^We^(*)|T(f))] = 


where |4/'(t)) = |\k(f)) and 


H'it) - i 


d 


|T'(f)) = 0, 


H\t) = e-^W. 

dt 

Here, we give the time-dependent generator S{t) of the transformation 
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(3) 
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with parameter ^ G [0,1] to be determined later 
Hamiltonian 


fl. 


Thus, we obtain the transformed 


H'(t) = -ujo i cos 
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—f sinfcnt) 
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(Tz + sin 


—f smlut) 
ui 


cr,, 


+—(1 - 0 cos{ujt)a^. 
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We make use of the identity exp sin(a;f)] = Yl°^=-oo (^0 exp(ma;f), in which Jn{-) is 
the nth-order Bessel function of the first kind, and divide the Hamiltonian into three parts 


Hq = ^UqJo ( ) cr^. 


H[{t) = ^(1 - 0 cos{ut)a^ + uqJi sin(a;f)cry. 


(7) 

( 8 ) 


= ojQ J2n ( — j j cos(2nn;t)(7, 


n=l 
oo 
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n=l 
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+ojo ^ J 2 n+i ( —^ 1 sin[(2n + l)ojt]ay, 
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where H^it) includes all higher-order harmonic terms (n > 2). Up till now, the treatment 
is exact. To proceed, we neglect the part H^lt) since the higher-order harmonic terms with 
the hi^er-order Bessel functions are reasonably negligible over a wide range of parameters 
space j^. Therefore, we arrive at the Hamiltonian H'(t) ~ i^Q -|- We verify that this 

treatment is beyond perturbation by the correct prediction of the BS shift for H/wq G [0, cxo). 

To proceed, we determine ^ by 


A 


OJ 


A 


A 


-? K'o = ^(1 -0 = -r- 


( 10 ) 


Consequently, we rewrite the Hamiltonian ff'(t) as 


H'{t) = 


-dwoff.+ 2(e ‘“V+ + e“‘(T_). 
oj J 4 


fl. 


( 11 ) 

Note that 


This is the counter-rotating hybridized rotating wave (CHRW) Hamiltonian 
the CHRW Hamiltonian possesses a RWA-like form with a renormalized transition frequency 
Jo OiJq and a renormalized driving strength A. The renormalized quantities in the trans¬ 
formed Hamiltonian results explicitly from the effects of CR interactions 
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We apply a rotating operation R{t) = exp (|a;tcr^) to the Hamiltonian (ITT]) , and get 


~ A H 


( 12 ) 


where A = Jq uq — ui is the effective detuning. The Hamiltonian H can be readily 
diagonalized. Its eigenenergies are 

(13) 

with hpi = ^A^ + H2/4 being the effective Rabi frequency. The corresponding eigenstates, 
i.e. dressed states are 


with 9 = arctan 


|±) = sin6'|=F) ± cos^li), (14) 

(flu — A)/(H/2) and |±) being the bare bases: cXzli) = ±|±). 


III. BLOCH-SIEGERT SHIFT 

A. Evaluation of the Bloch-Siegert shift 

We calculate the BS shift by using the effective Rabi frequency obtained in the former 
section. The BS shift measures the deviation of resonance frequency oj^es of the original 
Hamiltonian from that of the RWA; 


— Wes — W) 


(15) 


where ujq is the resonance frequency of the RWA case. In general, the resonance fre¬ 
quency cjres can be determined from the time-averaged transition probability in Ref. Q, 
P = [l — 4:{dqa/duo)'^] /2, where qa is the quasienergy of the Floquet Hamiltonian. In 
principal, on solving the equation dq^/duQ = 0 for driving frequency u, we determine the 
position of the resonant transitions. Alternatively, we can also use dq^/duQ = 0 as the 
resonance condition in the analytical calculation instead of dqa/duo = 0. 

Similar to Shirley’s method, the resonance condition in our formalism is dflj^/duo = 0 
since qa is related to Qr by the relation q± = (^u ± flR)/2. By solving the equation for 
variable u, we get the resonance frequency cUres- We can readily derive the explicit form for 
the derivative dflj^/duQ, which reads 


dnl 

dujQ 


OOqJo ( ^ 
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duo 


2x42(1-0 




duo' 


(16) 


where / Ouq is determined from Eq. ([TU]) and takes the form 

9? _ 2u,J. (4?) 


(17) 


9wo 7l {w + Wo [</o (fl;) - <72 (1^)]} 

Using Eqs. (fTOj) . (IT^ and (HU, we can self-consistently determine ^ and oj. The obtained 
u is actually the desired resonance frequency (Ures of fhe Rabi Hamiltonian. Therefore, it is 
straightforward to obtain the BS shift by Eq. ffTSj) . On one hand, in the weak driving case, 
we have correctly given the power series expansion of the BS shift in the previous work j^. 
On the other hand, it is easy to solve a large BS shift in the strong-driving limit. When 
H/cJo—)■! and d^/dojQ —)■ 0. Thus, the solution to Eq. flTBD for oj is determined 
by Jo (^0 = 0- Therefore, the hrst zero of Jq gives the resonance frequency in the 
strong-driving limit, 

Wres = x4/2.404826. (18) 

This result is the same as that obtained in Ref. Q. 

An interesting phenomena, coherent destruction of tunneling (CDT), occurs at A/u = 
2.404826 5|, l26|. Surprisingly, we immediately recognize that the CDT condition coincides 
with the resonance condition Eq. lfTSl) . Hence, it means that the CDT is actually a resonantly- 
driven phenomenon. In other words, it indicates that the large BS shift plays important role 
in this unique phenomenon. 


B. Numerical validation and comparison 

In order to illustrate the validity of our method and the accuracy of the calculated BS 
shift, we compare our results with the numerically exact and previous analytical results 
given by the methods based on Floquet theory. We show the approach of calculating the 
exact numerical BS shift by Floquet formalism. As shown in Ref. y, the derivative can 
be evaluated from the eigenstates of the Floquet Hamiltonian 

dqa 




^a^|( 7 /|A, 

■yl 


aO/ 


(19) 


where the notations are the same as Shirley’s (the index 7 represents states of the TLS 
and I denotes integer representing Fourier component. |Aqo) is the eigenstate of Floquet 
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Hamiltonian associated with quasienergy The coefficient = — (+)| if 7 represents the 
ground state (excited state)). For a fixed A and given u, we first numerically diagonalize 
the Floquet Hamiltonian (9) in Ref. Q with an appropriate truncation, and then get the 
eigenstate |Aco)- Using Eq. (IT^ we can numerically obtain the value of derivative 
with the fixed A and uj. It is therefore possible to determine the resonance frequency by 
numerically searching the solution to the equation = 0 or the minimum position of 
(^) • shift obtained in this way is referred to numerically exact result. 

We introduce two a ppr oximate methods for comparison by the analytical result of the q 
obtained by Shirley 2 , 19|. First, in terms of the analytical result ^ = 0, we obtain 

ojA? {2uo — oj)A^ 


CJ = CUo + 


4(a; + a;o)^ 64(a; + a;o)^ 

(9u;5 - 126cu^a;o + + A2uj‘^ujI - 2?>uuj^ - 8a;^)H® 


( 20 ) 


256(a; + Uo)^{9u‘^ — Uq)'^ 

and calculate the BS shift for comparison. On one hand, we can iteratively solve Eq. (|20|) for 
the resonance frequency u and gets the BS shift. It has been pointed out that this equation 
can be used to calculate the whole range of the shift from 0 to 100 % Q . On the other 
hand, carrying out perturbation calculation, we can derive the BS shift up to 6 th order in 
A, which reads 

M/zLU M/zLH 9.^(AIAA 

( 21 ) 


( 6 th) ^ (W ^ {A/ay 35(H/4)^ 


Wo 


4^0 


32ul 


In what follows, we show the comparison among various results. 

In Figs. [I](a)-[T](b), we make a quantitative comparison between the values of the BS shift 
obtained from our method and those given by the other approaches mentioned above over 
a wide driving strength range. From weak-driving to moderately strong-driving regime, we 
notice that our result is in good agreement with the numerically exact result and that of 
Eq. fl20|) . At the same time, it is obvious to see that the results of both 4th-order and 6 th- 
order BS shifts differ from those of the other methods when A/wq > 2. It means that the BS 
shift up to 6 th order of A is not adequate for moderately strong driving case and higher-order 
terms need to be considered when A increases. In addition, in Fig. [I](b), we find that our 
result is almost the same as the numerically exact result when A/wq > 5 but slightly differs 
from Eq. (|20|) when the driving is sufficiently strong [see the inset in Fig. [11(b)]. Moreover, 
it is obvious to see that the results given by Eq. (fT 8 |) provides a better description of the 
position of the resonance than that of Eq. (|20|) because the result given by Eq. (ITSl) becomes 
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more and more accurate as A —>■ cxo. 

In order to show the accuracy of our method in detail, we illustrate the quantitative 
difference of the BS shift among analytical and numerical approaches. First, we show the 
values of the BS shift calculated by the four methods in Table |T] for comparison. Second, 
in Fig|21 we show the deviation l^cugg — |in which denotes the 

numerically exact BS shift and represents the BS shift obtained by one of the three 
analytical methods, i.e., our method, Eqs. fl20l) and (flSj) . We hnd that when 2.5 < A/ojq < 
4.5, the deviation of our result from the numerically exact result is less than 1.2%, whereas 
the deviation of Eq. fl20p is less than 1%. Nevertheless, when A < 2.5a;o or A > A.duo, the 
deviation of our result is less than 1%. In particular, when Ajuj^ > 9, our method becomes 
the most accurate among the three analytical methods, with the deviation being less than 
0.02%. In addition, Eq. (ITSl) becomes more accurate than the Eq. (I20l) only for Aju^ > 17. 

To sum up, the BS shift given by our method is in quantitatively good agreement with 
the numerical result over the entire range 0 < Ajuj^ < cxo. Thus, it turns out that our 
calculation for the BS shift is beyond perturbation and the higher-order corrections involved 
in the BS shift have been correctly taken into account in our method. Moreover, our method 
is applicable to solve the dynamics of the TLS in the presence of the dissipation. As a 
result, we could demonstrate the measurable signatures of the BS shift in the excited-state 
population and the probe-pump spectrum. 


IV. SIGNATURES OF THE BLOCH-SIEGERT SHIFT 


In general, the signature of the BS shift can be measured in either the emission or absorp¬ 
tion process. In order to correctly discuss the observable signatures indicating the resonance, 
we should take into account the dissipation induced by the coupling to the environment in 
the Rabi model. We consider that the radiation of the TLS is described by the following 
master equation 

+p{t)a+a- - 2a-p{t)a+\, (22) 


where p{t) is the reduced density matrix of the TLS and k is the decay rate. 


the master equation has been obtained as a continued fraction in Ref. 


18|. 


The solution of 
Here, we apply 


8 



the unitary transformation to the master equation so as to transform it into a feasible form. 

By employing the relation between the transformed density matrix p{t) and the original 
density matrix p{t)\ 

p{t) = 


we derive the equation of motion for the element in the basis of |±), Pap{t) = (Q;|p(^)|/d): 

(23) 


d 


dt 

where a Greek index denotes ± and 

K 






n, A 


(24) 


with n taking on integer. The derivation of Eq. fl2^ is given in Appendix and the explicit 
expressions for ^ can be found therein. Moreover, the solution to the equation used 
in the following are given in Appendix [Bl Interestingly, we find that the derived master 
equation takes the similar structure as the Floquet-Born-Markov master equation given in 
Ref. j^l, which is derived from the total Hamiltonian consisting of the TLS, the reservoir, 
and their interaction. In the following, we calculate time-averaged population of TLS and 
the probe-pump spectrum based on the time-independent master equation we derived so as 
to show the detectable signals caused by the BS shift. 


A. The population of excited state of TLS 


Since the emission process is related to the steady-state population of excited state of the 
TLS, we attempt to find the signature of the BS shift by monitoring the population 18|. 


Moreover, the total intensity of resonance fluorescence is proportional to this population in 


the steady state 


In terms of the dressed-state population difference {Sz)ss = P-i-i-(oo) — 


p _(oo), we give the expression for the steady-state population of TLS by p(f). 


P++(oo) = lim Tr[|+)(+|p(f)] 

t^OO 


= lim Tr[i?(f)e'^*'*^-|-)(-l-e“ 

t^OO 

®“'fl*(<)p(i)l 

= ^ + ^(Sz)ss |cos(26') cos 

'A , ■ 

sin(a;f) 

OJ 
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+ sin(26') sin(ci;t) sin 
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^ sin(a;t) 


1 + (Sz}s 


cos(20)Jo ( ) + sin(26*) Ji ( — 
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OO 


n=l 


UJ 


2cOs{2e)J2n ( 


OJ 


sin(20) 




2n+l 


A 


UJ 


i — J2n-1 


A 


UJ 


cos{2nut). 


(25) 


When deriving the expression above, we have neglected the terms proportional to (s±)ss, 
which are of order k/Ci^. We hnd that the even multiple frequencies appear resulting from 
the unitary transformation. This property of the diagonal elements is similar as that of the 
ansatz used in Ref. 18| . 


Recalling that normal measurement of the experiment takes time average, we concentrate 
on the interested quantity, the time-averaged population /ujUjl from p++(cxd), which reads 


P++ = 


1 + {Sz)s 

jL 

2 k 7 ^’ 


A 


UJ 


cos(26')Jo -I-sin(26*) Ji 


A 


UJ 


(26) 


where we used 7o ^ k [cos(26') Jq sin(26') Ji (^^)] and (s^)ss ^ -lohz- 

In Fig. |3l we show the time-averaged population /ujTjT as a function of driving frequency 
UJ for various driving strengths A. The results show that the resonance is not only shifted 
but also broadened as A increases. These results are consistent with those of the previous 
work based on a continued fraction technique 18|. However, we point out that is 
physically related to the dressed-state population difference in our formalism, which is more 
clear than that in Ref. {^. More importantly, Eq. (|26|) indicates that the BS shift influences 
the dressed-state population difference. It means that probe-pump spectra may exhibit the 
effects of BS shift. In addition, according to these results, it is reasonable to predict that for 
the non-RWA case, the maximum intensity of resonance fluorescence can be achieved only 
when u = ujq + for a fixed A. 


B. The probe-pump spectrum 

In this section, we show the signature of BS shift from the probe-pump spectrum. It 
is well-known that the probe-pump spectrum is sensitive to the dressed-state population 
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difference. A non-zero population difference of dressed states can cause the lineshape of the 


spectrum to be asymmetric, which has two unequal sidebands 


29 


30| . In consequence of 


this property of the spectrum, we can expect that the spectrum may be sensitive to the BS 
shift because the BS shift can modify the dressed-state population. According to the linear 
response theory, the steady-state response of driven TLS to a weak probe held 

in which Qp and v are the amplitude and frequency of the probe held, is given by the Fourier 
transform of the two-time commutator 


3Q| 


S{v) oc Re lim / {[o'-it + r), cr+(t)]) 


t^OO 


(27) 


where o'±{t) = {t)a±U{t) are the operators in the Heisenberg picture with U{t) being 

the evolution operator for the total system consisting of the TLS, driving held and reser¬ 
voir causing the relaxation (in the absence of the probe held). We now apply our unitary 
transformation to the two-time commutator so that it may be related to the solution to 
Eq. (I23D. 

We hrst consider how to evaluate the two-time correlation function (cr_(t -|- r)(j_|_(t)) by 
making use of the unitary transformation and rotation operation because the other function 
((T_|_ (t)cT_(t -|- r)) can be treated similarly. The two-time correlation function is evaluated in 


the original frame as follows 


30|: 


{a_{t + T)a+{t)) = TTsR[U\t + T)a_U{t + T)U\t)a+U{t)p{0)pR] 

= TTsR[U{t)U\t + T)a_U{t + T)U\t)a+p{t)pR], (28) 

where p(0) is the initial state of the TLS and pr is the state of the reservoir. By employing 
U{t) = R(t + r)e®(*+")f/(t + r)17t(t)e-^WRl(t) (we assume that the evolution is r-dependent 
in the transformed frame) and with the aid of identities flAlSp . we have 


{a_{t + T)a+{t)) = Tr5R[i/t(r)R(t + r)e^(*+"V_e-^(*+")Rl'(t + r) 
xU{T)R{t)e^^*'> R\t)p{t)pR] 


n,l odd 


xf/(r)(Qe*'-* + Ae-*'‘"*)p(t)pH], 


(29) 
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where the summation is taken over all positive odd integers, and Cn and are defined in 
Eqs. (IA16j) and (IA17j) . 

To proceed, we take the long-time limit (t oo) and neglect the t-dependent terms in 
Eq. (|29|) because their contributions are negligible to a long-time observation. Besides, one 
can verify that the terms with the factor contribute a finite value to the spectrum only 
when u < 0. It is therefore reasonable to omit these terms when one considers the probe- 
pump spectrum. Consequently, we arrive at the following form of the two-time correlation 
function 


lim {a_{t + T)a+{t)) 

t^OO 


j Y, TrsR[UHT)ClU(r)C„f,(x,)p„]e 


nodd 


- Y 'IVs{ClTlK|C/(r)C„^(oo)p„C/t(r)]}e- 


nodd 




\p{0) = Cnp{oo) 


nodd 




(30) 


nodd 


where p(r) is the solution of Eq. (123|) with initial condition p(0) = Cnp{oo). When deriving 
the third line in Eq. flHOl) . we used quantum regression theory {mI- Similarly, the other 
two-time correlation function takes the same form as the last line in Eq. fl30|) but with initial 
condition p(0) = p(oo)(7„. 

Finally, we find that the two-time commutator is connected with the quantities in the 
transformed frame as follows: 


lim ([a_(f+ r),a+(f)]) 

t^OO 


3 ^ [(^l('^))lp( 0 ) = C'„p(oo) (^i('^))lp( 0 ) = p(oo)Cn 


nodd 


7 E«C»(^)))l«0).c.^(=o)-P(=o)c.. 


(31) 


nodd 


where ((Cl(r))) denotes mean value of averaged over the solution to the homogeneous 
part of master equation fl23|) with initial condition p(0) = C„p(oo) — p{oo)Cn- Substituting 
this expression into Eq. (HT} and letting p = —i{v — nu), we can express the spectrum 
function in terms of the Laplace transform of ((Cl(r))). Since ((Cl(r))) = /j/^,j((s_(r))) -|- 
f-,n{{^+{^))) + the spectrum function can be expressed as follows: 

^(u) oc ^Re [f+^g_{p) + /+„^+(p) + /+„c/(p)] \p=-i(v-no.), (32) 

nodd 
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where Qjijp) (j — ±, 2) are given in Eqs. (lR7[' - (lRr)[l and their initial conditions are given by 


xo ^ Tr{[s+,C„]p(oo)} 



(33) 


= -/+,v(Sr)ss + 

20 ^ Tr{[s;,C'„])i(oo)} 


(34) 


= 2fl„(i+U - 2/+„(S_) 


(36) 


Now we illustrate the role of BS shift in the probe-pump spectrum calculated without 
the RWA. In Figs. ID^a) andHU^c), we show the comparison between the non-RWA and RWA 
spectra for a hxed driving strength under the resonant condition u = ujq of the RWA. We 
notice that the RWA spectra are always symmetric with respect to the line oj = ojq while 
the non-RWA spectra are asymmetric and distinguish clearly from the RWA cases when 
A increases. In this case, the lineshapes for the non-RWA resemble those of the RWA at 
a detuning case [see Figs. Ill[b) andjU^d)]. In fact, this asymmetric behavior results from 
the BS shift. In Figs. Ill[b) andjlfd), we show that symmetric non-RWA spectra appear 
for u = Uq + 6ubs = ^res comparison with the RWA asymmetric lineshapes. These 
results indicate that when the BS shift is correctly compensated for driving frequency, the 
probe-pump spectra become symmetric and is similar as the RWA ones for u = ujq. 

In Fig. |5l we show the non-RWA spectra for A = O.lwo with three different pump frequen¬ 
cies. Although the detunings of the pump held are very small, one can observe difference in 
the lineshapes of the spectra. We notice that for uj = Wres = + (^i^bs = 0.000625ci;o), 

the curve exhibits two equally symmetric sidebands in the wing. However, for u 7^ (Ures, 
the lineshapes are clearly asymmetric (the spectra of a; = u^es ± <^i^bs are shown in the 
red-dashed line and the blue dotted-dashed line, respectively). Therefore, the probe-pump 
spectrum provides a way to sense the roles of the BS shift by two properties: (i) for u = Wq, 
the lineshape of the spectrum changes from nearly symmetric to asymmetric as A increases, 
while that of the RWA are always symmetric; (ii) the non-RWA lineshape becomes symmet¬ 
ric only under the exact-resonance condition, i.e. u = uJo + ^ojbs- One expects that the two 
properties can be checked in superconducting-circuit qubits under the moderately strong 
driving held. 
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V. CONCLUSION 


In summary, we have calculated the BS shift over the entire driving-strength range us¬ 
ing the derivative of the effective Rabi frequency and demonstrated detectable signatures 
induced by the BS shift. It turned out that our method can be correctly applied to study 
the small shift case when the driving strength is moderately weak, but also accurately solve 
the large shift case when the driving strength is sufficiently strong, which is beyond the 
perturbation theory. Moreover, our method allows us to examine the role of the BS shift in 
the emission and absorption processes from the open driven TLS. We showed that it is easy 
to obtain the time evolution of excited-state population of the TLS as well as its steady-state 
behavior, which indicates the emergence of resonance. We illustrated that the time-averaged 
population of TLS provides a direct measurement of the BS shift, which is consistent with 
previous work. Furthermore, we found that in experiment accessible parameters regimes, 
the non-RWA spectrum becomes symmetric only for u = uq + (Jcubs; he. for the resonance 
condition correctly taken into account the BS shift. Otherwise, the lineshape of the spectrum 
is generally asymmetric. While for ca = cuo, the non-RWA spectrum becomes asymmetric 
with the increase of the driving strength, while the RWA spectrum is always symmetric. 
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Appendix A: The derivation of master equation in the transformed frame 

We show the derivation of the master equation [Eq. (1231) ]. After the unitary transforma¬ 
tion and rotating transformation, we obtain the transformed equation as follows, 

+p(t)5-+(f)d_(t) -25-_(t)p(f)a-+(t)], (Al) 


14 


where o'±{t) = a±e Employing the eigenstates of H as the basis, we can 

obtain the equation of the elements of the density matrix 

d 

£(3)pa(3{t) - -{a\[a+{t)a_{t)p{t) 

+p{t)a+{t)a-{t) - 2a-{t)p{t)a+{t)]\P). (A2) 

To obtain the explicit form of the dissipation part of the equation, we insert the identity 
matrix such as l^)(^l = 1 between the operators, which leads to 

(d|d+(t)d_(t)p(t)|/3) = ^(d|d+(t)|i>)(i>|d_(t)|/f)(/f|p(t)|/3) 

= (A3) 

n,n' 

fi,!/ n,n' 

^l,v n,n' 

where are the Fourier coefficients from the expansions {pi\a±{t)\(5) = 

To proceed, we need to give the explicit forms for which can be evaluated by the 

following integral 

^a/3,n = ^ dt{ua{t)\a±\up{t))e-^^‘^^. (A6) 

It is straightforward to calculate the expressions for which are given by 

xp.n = 51 ](hA«+ (A7) 

I odd 

x:-.„ = 5 + f:A-A (as) 

I odd 

A'-+,„ = lY.(f-.An + f;A,-nl (A9) 

I odd 

= -i Y,(PA- + A.A-A (Aio) 

I odd 

where the summation is taken over all positive odd integers, and 

f+^i = - ± A _1 cos^ d =F Ji+I sin(26'), (All) 

sin(26*), (A12) 

fti = ^ ± di-i T Ji+i sin(20) T Ji cos(20). (A13) 


f^i = 5i^i ± A_i —i sin^ 0±Ji 
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The explicit expressions for ^ can be directly obtained by the relation 

(AM) 

When calculating we used the identities 

= \Y^ (A15) 

n odd 

where 


Cn — /+,n'5+ + (A16) 

Dn = fZ,nS+ + /+,nS- + f~nSz- (Al7) 

Here, we introduced a set of dressed-state operators: 

s+= sL = |+)(-|, S 2 : = |+)(+| - (A18) 


It turns out that the master equation in the new representation still possesses explicit 
time-dependence. However, in the strong-driving regime, it is feasible to remove the explicit 
time-dependence by invoking the partial secular approximation (Moderate RWA) 3|, l27|. In 
other words, we only keep the terms satisfying n -|- n' = 0 in Eqs. (lASP - flASp . This approx¬ 
imation can be justified when Qr 3> k. Consequently, we arrive at the time-independent 
master equation in Eq. fl23p . One can verify that the time-independent master equation 
can predict almost the same dynamics as that given by Eq. (|22l) when Qr 3> k. 


Appendix B: The solutions to the master equation in the transformed frame 

The master equation can be rewritten in terms of the mean value of dressed-state oper¬ 
ators by the following relations; 

= p++(t) - p—{t), {s+{t)) = {s-{t))* = p-+{t), (Bl) 

I {i.m = -7. {S.(*)) - 27.((i+(()) + {i-m - 70 , (B2) 

^ jps-(t)r 

= iVtR (s+(t)) - 71 {s,{t)) - 7_ (s_(t)) - 7+ (s+(t)) - 72, (B3) 
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where 


Iz 

= '^++,++ 

“ '^++,-) 

7o 

= '^++,++ 

+ ^ + + , -5 

7i 

= '^++,-1—: 

) 

72 

— ('^--1-,++ + ^- + , - 

7- 


> 

7+ 




(B4) 

On solving the equations, we can fully determine the evolution of the TLS. In particular, 
the steady-state solutions can be easily found as follows: 

- 0^70 - 47172(7- - 7+) + 7 o( 7 - - ll ) 


('^^)ss ~ 


(5+)ss = (s-)ss = 


47?(7-- 7+) + - 7-+ 7i)72 

{i^R - 7- -h 7+)(7 o7i - l2lz) 


(B5) 

(B6) 


47 ?( 7 - - 7+) + - 7- + iDlz 

When calculating the probe-pump spectrum, we need the solutions to homogeneous parts 
of Eqs. (|23|) . The corresponding homogeneous differential equations can be solved by Laplace 
transform. Denoting the solutions as {{sj{t))) {j = ^,z), we can obtain their Laplace 
transforms, 

noo 

9+{p) = / e“P*((s+(f)))df 

Jo 

1 


{xoiip + 7+ + i^R){p + Iz) - 27i] + ?/o[27i - 7-(p + 7^)] 


F{p) 

-liZo{p + - 7- + 7+)}, 

poo 

9-i.p) = / e-P\{s-{t)))dt 

Jo 

1 


(B7) 


F{p) + 7z) - 27 i] + 2:0 [27i - 7-(p + 7^)] 

-liZo{p - i^R - 7- + 7+)}, 

poo 

9 zip) = / e-P^{{s,it)))dt 

Jo 

= 7^{^o[(p + 7+)^ + S2 r - 1-] - 27iPo(p + - 7- + 7+) 

-27iPo(p - - 7- + 7+)}. 


(B8) 


(B9) 


where the initial conditions are Xq = Tr[s+p(0)], Hq = Tr[s_p(0)] and Zq = Tr[s^p(0)], and 
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the polynomial F{p) is given by 


F{p) = p^ + 472(7- - 7+) + (fi| - 7 - + iDlz 

+P^{lz + 27 +) + p{h\ - 472 - 72 + ^2 ^ 27 + 7 ^). (BIO) 
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AjujQ numerical present paper Eq. (I2UI) Eq. (I18p 


1.0 

0.063224 

0.063268 

0.063228 


3.5 

0.707959 

0.716200 

0.712320 

0.455407 

6.0 

1.641809 

1.649924 

1.650482 

1.494983 

8.5 

2.637787 

2.640075 

2.639255 

2.534559 

11.0 

3.653740 

3.652351 

3.641373 

3.574136 

13.5 

4.678502 

4.675271 

4.650384 

4.613712 

16.0 

5.707919 

5.703825 

5.664602 

5.653289 

18.5 

6.740093 

6.735637 

6.683190 

6.692864 

21.0 

7.774035 

7.769474 

7.705492 

7.732441 


Table I: The comparison of the BS shift obtained by various methods. 



Eigure 2: (Color online) The deviation between analytical results and numerical result. 


20 











Figure 3: (Color online) The time-averaged population as a function of driving frequency u 
for K = 2 X 10“^a;o and various driving strength A. 
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Figure 4: (Color online) The probe-pump spectrum S{v) is shown as a function of probe frequency 
r; for K = 2 X 10“^a;o- 
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Figure 5: (Color online) The probe-pump spectrum S(v} is shown as a function of probe frequency 
V for K = 2 X 10“^a;o with A = O.lwo- 
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